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1. Introduction: Wilson twisted mass action 

In this contribution, we study the cutoff effects of observables computed on a lattice with lattice 
spacing a when Wilson twisted mass fermions at maximal twist are considered 1 . The twisted mass 
QCD action in the continuum, at tree-level of perturbation theory, is given by 

Sf\X,X} = j d 4 xx{x){Y^+m + in q Y 5 T 3 )x(x), (1.1) 

with T 3 the third Pauli matrix acting in flavour space and {%,%} the so-called twisted basis. In 
equation ( |1.1[ ), mo denotes the bare untwisted quark mass and \x q the bare twisted quark mass. The 
mass term can be written in terms of a polar mass M and a polar angle a as 



m Q + in q y 5 z 3 = M e iay5t3 with M= Jml + n* a = arctan(^). (1.2) 

V H mo 

The twisted basis {%,%} is related to the so-called physical basis {y, ij/} by the non-anomalous 
axial transformation 

^( x ) = e ia> y^x{x) xj/(x) = x{x)e ia ^ , (1.3) 

for the particular choice of the twisting angle (0 = a, since it brings the twisted mass QCD action 
back to the standard form. The Wilson-regularized twisted mass action (Wtm) written in the twisted 
basis has the form 

ST m lX,X] = « 4 I X(x) [D w +m + i W z 3 ]x(x) ; D w = hy^ + d^-ard^} (1.4) 

x Z 

with Dw the Wilson-Dirac operator of the free case, dp and d* the forward and backward partial 
lattice derivatives and — y^u^h tne Wilson term. 

2. Wilson twisted mass free-fermion propagator with infinite time-extent lattices 

The expression for the Wilson twisted mass fermion propagator in the twisted basis, at tree-level 
of perturbation theory (PT) and in momentum space is given by 

~ = -/Isin( a p4)/l/-^l/+[^(l-cos(ap 4 ))+^]ll / -/^y 5 ^ (2 {) 
^sm 2 (ap 4 ) + ^ 2 + ^(l-cos(a/? 4 ))+-#] 2 + jUf 

where 1 and 1 f are the identity matrices in Dirac and flavour space. The structure in colour space 
has not been written since it is just an identity matrix at tree level of PT and we have defined, 

J({p) =m + -f^ sin 2 (^i), J(r=lY r i sm(ap i ). (2.2) 
a £ 2 a ( tt 

We obtain the expression for the quark propagator in the time-momentum representation when an 
infinite extension of the time direction is considered. To perform the integral in the continuous 



The interested reader may find a much detailed account of this work in the thesis of J. G. L., see ref. 



2 



Cutoff effects for Wilson twisted mass fermions at tree-level of perturbation theory 



J. Gonzalez Lopez 



momentum P4 amounts to performing a contour integral and compute the residues of the integrand. 
The integration contour encloses the poles of the integrand at the energy points 

( p 4 = \iE x cosh£i = ^~Ky*) physical pole 

\ P4 = 1 -(x + iE 2 ) cosh E 2 = doubler pole 

with ^ 2 = (1 +arJt) 2 + (l - r 2 ) (a 2 J^ 2 + a 2 ^ 2 + a 2 ^). The final expression for the propa- 
gator, in the limit r — ► 1, is then 2 , 



S„(p,t) = 2 sm h e { s § n (0 sinh g i 74 V ~ V + [(1 - cosh£i) 11/ 



ia^}e- E ^\ + 8 ifi 1 11 



(2.3) 



where coshfij = 1 + 2(1W) 



The fermion propagator is a matrix in Dirac space and can hence be decomposed in terms of the 
Dirac gamma matrices as Soo(p,t) = Su(p,t) 1 + £ M S^(p,t) Y^+S 5 (p,t) y 5 . 

3. Hadron correlation functions 

3.1 Pseudo-scalar meson 

The interpolating fields describing the charged pions, n + and n~ respectively, in the physical and 
the twisted bases are 

^(x) = &>\x)^i& >2 (x)=P\x)^iP 2 (x) (3.1) 

where £P a (x) = \\f(x)y$^\\f{x), with a = 1,2,3, is the pseudo-scalar density written in the physical 
basis while P a (x) = X(x)y^x{ x ) i s the pseudo-scalar density written in the twisted basis. 
The time dependence of the two-point correlation function for the charged pseudo-scalar meson in 
the time-momentum representation is given by 

C^ { t) = ^ J E{\S" u (p,t)\ 2 +i \S u ,(p,t)\ 2 + \S u 5 (p,t)\ 2 }. (3.2) 
P M=l 

We denote the Wilson twisted mass fermion propagator in the twisted basis for "w quarks" as 5". 
Note that this is not identical in the twisted mass case with the propagator for "d quarks" denoted 
as S d . N c (Nd) is the number of colours (Dirac components). 

3.2 Proton 

The local interpolating field describing the proton in both bases is given by 3 , 

&> a (x) = -V2e ahc [d^(x)C- 1 Y 5 u h (x)} u a Ax) = V2e ahc [d T a {x)C- { Y 5 u h {x)] e^Ua^x). (3.3) 

2 This expression is valid for all possible values of the discrete Euclidean time t, including if it is negative or zero. 
The function sgn(f ) is the sign of t, and we have denoted sgn(0) = 0. It is just a convention in order to give one general 
expression for the propagator for all possible values of t. 

3 The Greek (Latin) letters denote Dirac (colour) components and u, d denote the flavour content. The notation used 
for the flavour structure is y = (j) and %= QJ) . C is the charge conjugation matrix and [] denotes spin trace. 



3 



Cutoff effects for Wilson twisted mass fermions at tree-level of perturbation theory 



J. Gonzalez Lopez 



The expression for the time dependence of the proton correlation function is then 

C^ff, (/) = -jj^- ££{cos(a>) L v (p, q,t)+L 4 (p,q,t) + isia(co)L 5 (p, q,t)}, (3.4) 
p q 

with the definitions 

Lu,5 (p,q,t)=S u Uj5 (-(p + q),t)\(N d + l)S u u (p,t)S l {j (q, t) + (N d + 3) £ (p,t)S» (q, t) 

1 n=i ' (3.5) 

-(N d + l)S$(p,t)S$tf,t)} 

L»(j},q,t)=Sl(-(j}+q),t)\ (N d + 3)Sf,(p,t)S u u (q,t) + (N d + l) £ S£(jM)S£(?,0 

1 m=i ' (3.6) 

-(N d + 3)S u 5 (p,t)S u 5 (q,t)}- 

4. Scaling test 

4.1 Wilson average (WA) and mass average (MA) for standard Wilson fermions 

In Ref. it has been demonstrated that when averaging physical observables computed with 
Wilson actions having opposite signs of the quark mass (MA) or opposite signs of the Wilson 
parameter (WA), these quantities are 0(a) improved. Since WA and MA are equivalent, we show in 
Figure |], only the cutoff effects at the example of the proton mass when the MA is performed. In the 
left graph the behaviour of the proton mass NM P as a function of ^ is given when a standard Wilson 
regularization is used. In order to describe the behaviour of the physical quantities computed 
analytically at selected values of 4, we use the following fitting functions: 

yi = a + a^ + a 2 ^ y 2 = b + + bl ^4 • ^ 

Here y\ (y 2 ) is the physical observable under consideration and its value in the continuum limit is 
given by the coefficient oq {bo). We use two functional forms, the first formula of equation ( fi~T| ) 
for a leading 4 behaviour (standard Wilson fermions) and the second formula for 0(a) -improved 
quantities. 



The two lines in the left graph of Figure [j] originate from a fit to equation (4. 1) and correspond to 
the proton mass obtained from the same Wilson actions differing only in the sign of the quark mass. 
The linear behaviour in 4 shows the 0(a) scaling violations present in the standard Wilson theory. 
From the plot it is clear that in both cases the value of the proton mass in the continuum limit 
is the same and the expected one at tree-level of PT. From the fit, the corresponding coefficients 
a\ turn out to be the same in magnitude but have opposite signs for positive and negative quark 
masses. Thus, performing the (MA), it is to be expected that the 0(a) effects cancel and the 
scaling behaviour changes drastically from a | to a p behaviour. This can indeed be seen in 
the right graph of Figure [j] Inspecting the fit coefficients a 2 and fej, we find a 2 w f« bf A as 0.5. 
Therefore, the magnitude of the leading order cutoff effects does not only change from an 0(a) 
to an 0(a 2 ) behaviour but also the 0(a 2 ) do not increase on performing the Wilson average with 
respect to the standard case. 
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Figure 1: In the left graph, the cutoff effects and the continuum limit of the proton mass obtained from two 
standard Wilson actions differing only in the sign of the quark mass, |iVffio| = 0.8 are shown. The lattices 
are 4 < N < 20. In the right graph, the average of the proton masses obtained from the same two standard 
Wilson regularizations with quark masses A^mo = ±0.8 (MA) has been calculated. 




1/N 2 

Figure 2: Cutoff effects in the proton mass with Wtm fermions at maximal twist. To study the case of 
maximal twist, we set the quark mass to N[l q = 0.5. The lattices considered are taken as 4 < < 24. 



4.2 Wilson twisted mass fermions at maximal twist 

Instead of performing a MA or WA, a way to obtain an automatic 0(a) improvement is to work 
with Wtm fermions at maximal twist Maximal twist is reached for a value of the twist angle of 
(0 = n/2. At tree-level of perturbation theory this can be achieved by simply setting the untwisted 
quark mass mo = 0. The value of the quark mass is now fully given by the twisted quark mass, 
NM = Nfl q . The proton mass as a function of is shown in Figure |2| and the corresponding fit is 



performed using the second fitting function of equation (|4JJ). The figure clearly demonstrates that 
automatic 0(a) improvement is indeed working and that the cutoff effects have changed from a jj 
behaviour of standard Wilson fermions to a ^ behaviour when maximally twisted mass Wilson 
fermions are employed. Moreover, as a result of the fit, the coefficient b\ comes out to be very 
small, bf = 0(l(r 2 ). This value is one order of magnitude smaller than the corresponding coef- 
ficient a2 of the 0(a 2 ) effects for standard Wilson fermions which we find to be = O(l0~ l ). 
Note the the value of b\ is also smaller than the one for the case of MA discussed above. 



5 



Cutoff effects for Wilson twisted mass fermions at tree-level of perturbation theory 



J. Gonzalez Lopez 



We have also analyzed the chiral limit behaviour of the pion and proton masses by determining the 
coefficients a\ and b\ of equations (fO|). In the case of standard Wilson fermions, the coefficient 
a\ which determines the size of the 0(d) cutoff effects vanishes in the chiral limit thus leading 
to only 0(a 2 ) scaling violations in the massless theory. For Wtm fermions at maximal twist the 
situation is even better. Here, the coefficient of the 0(a 2 ) cutoff effects vanishes in the chiral limit, 
thus leading to scaling violations of 0(a A ) only since all odd powers of a vanish for maximally 
twisted mass fermions. Therefore, for Wilson twisted mass fermions at maximal twist the breaking 
of chiral symmetry at finite lattice spacing is much smaller than for standard Wilson fermions and 
maximally twisted mass fermions are indeed chirally improved. 

4.2.1 Out of maximal twist 

Here we want to study a situation when we allow an 0(a) error in setting the untwisted quark mass 
to zero. In order to realize this situation at tree-level of PT we 'force' these effects by simply fixing 
the twisted mass to be the physical quark mass and the untwisted mass is set to be proportional to 
jf, as N/Aq = a and Nmo = jj ^ 0(a) where a is kept fixed and j8 is a measure parametrizing the 
amount of violation of the maximal twist setup. The twist angle CO and the bare polar mass M can 
be obtained as a function of a and /3 as 

Therefore, even if the condition of maximal twist can be only obtained up to 0(a) cutoff effects, 
which is generically the case in practical numerical simulations, the observables, which are only 
functions of the polar mass, are still automatically 0(a) improved. 

Moreover, equation ( fO| ) also shows how the size of the leading discretization effects depends on 
the ratio between the untwisted and twisted quark masses. This ratio in turn determines the value of 
the lattice spacing at which the asymptotic 4r scaling sets in. Only when this ratio is small enough 
and hence the lattice does not need to be chosen too large a reliable continuum limit using reason- 
ably sized lattices can be performed. The left graph of Figure || demonstrates that the asymptotic 
scaling sets in for lattices with 4 <N <64 when < ^ < 2. 

However, for ^ > 10 the continuum limit is not reliable anymore if N is chosen to be too small. 
This can be seen in the right graph of Figure [|. Using only small values of N leads to an inconsistent 
continuum limit value. Therefore, larger lattices are needed in order to obtain the correct continuum 
behaviour as can be also seen in the right graph of Figure [| Here we have added a fit of the data 
for a value of /3 = 10.0 taking only large lattices into account, i.e. using only values of N > 40. In 
this case, indeed the right continuum value is obtained. Of course, for practical simulations, using 
only lattices with N 3> 40 appears to be rather unrealistic. 

5. Conclusions 

In this contribution we have demonstrated at tree-level of perturbation theory that when Wtm 
fermions at maximal twist are considered, physically relevant quantities are automatically 0(a) 
improved. In addition, the magnitude of the leading 0(a 2 ) corrections is rather small. In the chiral 
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Figure 3: Left graph: Behaviour of the pion mass as a function of jL, for lattices with size 4- <N < 64. The 
twisted quark mass is set to Njl q = 1.0 and the untwisted quark mass is zero up to 0(a) cutoff effects i.e. 
Nmo = jj with p = 0.0, 1.0,2.0, 10.0. Right graph: a zoom of the graph on the left with an additional fit for 
the analytical data corresponding to j3 = 10.0 which considers only large lattices 40 <N < 64. 

limit the 0(a 2 ) (0(a)) effects disappear leading then to scaling violations of 0(a 4 ) (0(a 2 )) in the 
case of Wtm fermions at maximal twist (Wilson fermions). Therefore, Wtm fermions at maxi- 
mal twist show a substantially improved scaling and chiral behaviour when compared to standard 
Wilson fermions which render Wtm fermions a powerful formulation of lattice QCD. 
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